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Abstract 


The  paper  treats  the  general  theory  of  optimal  confidence 
limits  for  systems  reliability  introduced  by  Buehler 
These  general  statements  are  specialized  to  the  case  of  series 
systems.  It  is  noted  that  many  results  previously  given  are 
false.  In  particular,  counterexamples  for  results  of  Sudakov^ 


$2 - H9T4>^  Winterbottom  (1974)  and  Harris  and  Soms  (1980,1981)  are 

given.  Numerical  examples  are  provided,  which  suggest  that 
despite  the  deficiencies  of  these  results,  they  are  nevertheless 
valid  for  those  significance  levels  likely  to  be  used  in 


practice. 
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1.  Introduction  and  Sunnary 

A  problen  of  substantial  inportance  to  practitioners  in 
reliability  is  the  statistical  estination  of  the  reliability  of  a 
systen  of  stochastically  independent  coaponents  using  experiaental 
data  collected  on  the  individual  coaponents.  In  the  situations 
discussed  in  this  paper,  the  conponent  data  consist  of  a  sequence 
of  Bernoulli  trials.  Thus,  for  conponent  i,  i«l,2,...,k,  the 
data  is  the  pair  (n^.Y^),  where  n^  is  the  nunber  of  trials  and  Y^ 
is  the  nunber  of  observations  for  which  the  conponent  functions. 
*1**2*‘‘* **k  are  *ssun«d  to  be  nutually  independent  randon 
variables. 

This  problen  was  treated  in  Sudakov  (1974),  hint erbot ton 
(1974),  and  Harris  and  Sons  (1980,1981);  one  purpose  of  the  pre- 
sent  paper  is  to  exhibit  coUnterexanples  to  theorens  in  the  above 
papers. 

In  Section  2  we  discuss  the  general  theory  of  optinal 
confidence  Units  for  systen  reliability  so  that  the  notation  and 
definitions  to  be  enployed  in  the  balance  of  the  paper  have  been 
prescribed.  Sone  general  results  on  optinal  confidence  linits 
are  established. 

In  Section  3  the  counterexanples  previously  nentioned  are 
exhibited  and  the  specific  errors  in  the  proofs  of  the  theorens 
are  indicated. 

Section  4  presents  the  proof  of  a  special  case  of  the  key 
test  theoren  (Ninterbotton  (1974)),  the  general  forn  of  which  was 
invalidated  by  a  counterexanple  in  Section  3. 

The  consequences  for  reliability  applications  are  discussed 


in  Section  S. 
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2.  Bushier' s  Method  for  Optical  Lower  Confidence 
Bounds  for  Systea  Reliability 

We  now  introduce  the  notation,  definitions,  and  assuaptions 
that  will  be  used  throughout  the  balance  of  this  paper. 

1.  Let  pi#  i»l,2,...,k  denote  the  probability  that  the  itk 

coaponent  functions.  The  coaponents  will  be  assuaed  to  be 
stochastically  independent.  The  reliability  of  the  systea  will 
be  denoted  by  h(p)  ,  where  p  ■  CPj*P2*  •  •  •  *Pjc)  *  0  <  <  1.  It  is 

assuaed  that  h(0,0,...,0)  *  0,  h(l,l,...,l)  ■  1,  and  that  h(p)  is 
non-decreasing  in  each  p^,  i«l,2,...,k.  Further,  h(p)  is 
continuous  on  (p|0  <  p^  <  1},  which  follows  readily  froa  the 
assuaption  of  independence.  These  properties  hold  for  coherent 
systeas  (see  Barlow  and  Proschan  (1975)). 

2.  Let  S  «  (x|xi-0,  1,...;  n^  i»l,2, . . , ,k).  g(x)  is  said  to 
be  an  ordering  function  if  for  Xj  <  z^,  x2  <  z2,...,  xk  -  xk* 

x,  i c S,  g(x)  >  g(z).  (It  is  often  convenient  to  noraalize  g(x) 
by  letting  g(0)  -  1  and  g(nj ,n2, . . . .n^)  *  0.  With  such  a 
noraalization,  g(x)  is  often  selected  to  be  a  point  estiaator 
of  h(p).) 

3.  Let  R  •  {rj,r2, . . . ,r#,  s  >  2>  be  the  range  set  of  g(x).  With 
no  loss  of  generality  we  order  R  so  that  r j  >  r2  >  . , .  >  r# . 

4.  Let  A^  ■  (xjg(x)  ■  r^,  x * S,  i-1,2, . . .  ,s).  The  sets 
constitute  a  partition  of  S  induced  by  g(x). 

5.  We  assuae  throughout  that  the  data  is  distributed  by 


k  fnil  Bi_xi  xi  x  fni)  rl  Bi"yi 

f(xjp)  -  Pi(X-x)  -  n  1  pii  1  -  n  jpt  R* 

P  i-l^x^  i«l 'y^'  * 


*  rni)  yi  "i-yi 


where  q4  ■  l-pi# 


n^-y^ ,  i*l , 2, . . . ,k. 


With  no  loss  of 


(2.1) 
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generality,  we  assume  nj  <  n2  <  . ..  <  nk. 

From  these  definitions,  it  follows  that 

pp{x  c  j,  Ai)  ‘  i  ri)  •  **•« 

From  (2.1)  and  (2.2),  we  have 

U1  u2  uk 

l  l  ...  I  *C*;p)  ,  (2.3) 

ij*0  ij^O  i k» 0 

where  i  ■  (ij ,ij , . . .  ,1^)  and  u2  •  u2 (i i)  » •  •  •  »ujt  *  “jt(i1»i2»  •  *  • » 
ij^j)  are  integers  determined  by  r j . 

6.  Subsequently  we  will  need  to  extend  the  definitions  of  S  and 
g(x)  to  real  values.  We  denote  this  as  follows.  Let 

s*  •'  {*1°  5  *i  5  V  1"1"2 . k} 

and  let  g'(x)  be  defined  on  S*  with  g(x)  >  g(z),  x,  i  e  S* ,  when- 
ever  x  <  z,  and  g’(x)  «  g(x)  for  xeS. 

Then 

.  ,  «*1>  It2]  Ctk] 

Pd1*CX)  >  r.[  -  l  l  ...  I  f(I;p)  ,  (2.4) 

Pl  **  ij»0  i2«0  ik»0 

where  t2  ■  t2(ij)  , . .  .  ,tk  ■  tj?  (ij  ,  i2 » •  •  • » *k_  j)  >  with 

tx  ■  sup|t j t  e  S*  and  g(t,0,0, . . .  ,0)  >  r^jand  tk(i1,i2» • . • .ij.j) 

*  supjt 1 1  e  S  and  g(ij,i2,..«,ik_^,t,0,..,0)  >  r^| (  t*2,3,...,kl 
We  now  introduce  the  notion  of  Buehler  optimal  confidence 
bounds.  Let  g(x)  ■  r . .  Then  define 

*g(x)  "  inf{hCP)iPp{JIg(i)  >  g(x)}  >  a}  .  (2.5) 

Equivalently,  by  (2.2),  we  can  also  write 


By  assumption  1,  h(p)  is  continuous  and  the  set  of  parameter 
points  satisfying  (2.5)  is  compact;  therefore  the  infiaua  in  (2.5) 
and  (2.6)  is  attained. 

Assume  that  there  is  an  integer  j,  1  <  j  <  s-1,  such  that 
br  >  ar  .  Then  there  exists  a  pQ  such  that 


>  a 


TJ 


i 

u 

i-1 


(2.7) 


In  Addition,  there  exists  e  such  that 

p5  (xe  u  A,}  >  o  ,  h(p.)  <  b  . 

pll  i-i  xi  1  rj 

Since  b  >  b  >  ...  >  b  ,  fron  (2.7)  we  have 
1  ‘  2  ‘  s 


(2.8) 


b(p.)  <  b  ,  1  ■  1, 2,  , . . ,  j  . 

1  rt 


Therefore 


(2.9) 


\{u  A  Ai}  s  pfl{ht5i)  <  b«tx)}'  (2-10> 


which  is  a  contradiction.  Consequently,  there  is  no  integer  j, 
i  <  j  <  s - 1 ,  for  which  b_  >  a  . 

'  *  .  Ti  Ti 

Remarks.  Fron  (2.6),  it  follows  that  a_  ■  0  and  b_  is  also 

-  rs  s 

necessarily  2ero.  Note  further  that  in  (2.7)  it  is  possible 

that  the  infinua  is  attained  at  a  point  for  which  P~|x  e  U  AjJ  > 


To  see  this  consider  the  following  example. 


Let  k  ■  2,  nj  ■  5,  n2  ■  10,000,  Xj  -  0,  x2  ■  5, 

g(x)  *  nj+nj-Xj-Xj,  h(p)  ■  PjP2*  Jt  is  ®*»il Y  seen  that  the 

hypotheses  of  Theorem  2.1  are  satisfied.  Thus,  for  the  data 

6 

given,  g(x)  ■  10,000  ■  r,.  The  set  U  A,  consists  of  all  points 

1-1  r  l 

(Xj,x2)  for  which  Xj  +  x2  -  s*  th,t  is*  A1  * 

A2  "  jc1*0) » (®» 1)|»  *nd  50  on*  Consequently, 


S’W*',!,  *4  - 


includes  the  parameter  points  (0,  P2a)  where  p2a  satisfies 

??  {*2*°}  -  a*  ,ince  Pp  {xi  ~  5}  *  1  when  P1  *  °*  TIlu5  in{ 

h(p)  -  0  for  all  0  <  a  <  1. 

The  reader  should  also  note  that  the  aonotonicity  of  h(p)  is 
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not  utilized  in  the  proof,  which  is  valid  whenever  h(p)  is 
continuous . 

It  is  easy  to  see  that  a.-,,  is  aonotone,  i.e.,  a„  >  a  > 

g(xj  rl  "  r2 

...  >  ar  .  This  follows  froa  (2.7)  upon  noting  that  as  j 
increases,  the  set  of  p  satisfying  (2.7)  increases  and  the 
infiaua  is  taken  over  a  larger  set. 


k 

Corollary.  For  a  series  systea  h(p)  -  II  p4 .  Then  if 
- k  k  i-1 

g(x)  »  H  (n. -x. )/n.  «  n  y,/n4,  the  hypotheses  of  Theorea  2.1 
i-1  11  1  i-1  1  1 

are  satisfied  and  the  conclusion  follows. 


k 

Reaark.  Mote  that  g(x)  -  U  (n4-x,)/n<  is  the  aaxiaua  likelihood 

i-i 

estiaator  as  well  as  the  ainiaua  variance  unbiased  estiaator  of 
k 

n  p.  and  is  therefore  a  natural  choice  of  an  ordering  function 
i-1  1 

for  this  case. 

We  now  establish  the  following  theorea. 


Theorea  2.2.  Let  g(x)  •  and  let 

f  (x;a)  -  sup  P-fg(X)  >  r.},  0<a<l  .  (2.10) 

h(p)-a 

Then 

sup  f*(x;a)  -  1 
0<a<l 

and  f*(x;a)  is  non-decreasing  in  a. 

Proof.  Since  h(p)  is  continuous  and  h(l)  -  1, 

lia  sup  P5{g(X)  >  r.|  ■  1  . 
a-*-l  h(p)-a 

Now  choose  a  and  b  such  that  0<a<b<l, 

i  ri}  ‘ 

and 
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>  r^J  •  f#(x;b)  . 

Let  be  the  set  of  indices  i  such  thet  p^a  <  1.  Then  it  is 
possible  to  replace  pia  by  p’b,  ie  Ia,  where  pla  <  p{b  <  1,  so 
that  htPjP  ■  b,  where  p^b  -  pia,  ie  I*.  This  follows  since 
h(l)  *  1  >  a  and  h(p)  is  continuous.  The  conclusion  follows  from 
the  aonotone  likelihood  ratio  property  of  the  blnoaial  distribution 

Remark.  Only  the  continuity  of  h(p)  was  used  in  the  proof  of 
Theorem  2.2. 

Por  the  case  of  series  systems,  it  is  possible  to  strengthen 
Theorem  2.2  and  to  exhibit  the  above  construction.  This  is  done 
below. 

k  # 

Corollary.  Let  g(x)  «  r  If  h(p)  ■  n  p  ,  then  inf  f  (x;a)  «0 
*  3  i*l  0<a<l 

and  f  (x;a)  is  strictly  increasing  in  a  whenever  all  u^  <  n^ 

(see  (2.3)  for  the  definition  of  u^),  j»l,2,...k. 

Proof.  From  the  hypotheses, 

i  rj}  i  '-’i"1  •  . k  • 

k 

and  since  n  p.  -*•  0  implies  at  least  one  p.  0,  this  gives 
i-1  1  1 

inf  f  (x;a)  »  0  . 

0<a<l 


*  - 

To  show  that  f  (x;a)  is  strictly  increasing  in  a,  consider 

0<a<b<l  and  let  pa  -  (pgl . pafc)  satisfy  f*(x;a)  »  P-  |g(X)  >  r^ 

Similarly,  let  pb  satisfy  f#(x;b)  ■  P~  jg(X)  >  r^J.  Let 
Ia  •  jij.ij,. . . ,ir|  be  any  non-empty  set  of  indices  such  that 


fb. 1/r 


p.i  (— )  <1  (non-empty  because  otherwise  multiplying  the 

*  J  * 


i 


,  J 
? 


components  would  give  b  >  1,  a  contradiction)  and  let  I®  be  the 
remaining  indices.  Then 

b  1/r 

C  n  P4i  (|)  )  n  c  Pai  -  *  •  (2.ii) 

j€la  alj  a  jel*  aij 

From  the  monotone'  likelihood  ratio  property  of  the  binomial 
distribution, 

p5.{*c«  :  rj}  <  Pp-(stx)  >  Tj}  . 

where  the  components  of  p*  are  given  by  (2.11).  This  gives 

f#(x;a)  <  f* (x;bD  , 
which  is  the  desired  conclusion. 

Remarks.  Note  that  if  at  least  one  u,  •  n.,  it  follows  immedi- 

-  3  J  k 

ately  from  C2.S)  that  hCp)  »  0.  For  g(x)  ■  II  (n,-x. )/n,  the 

i»l  1  1  1 

condition  u^  <  n^  is  equivalent  to  x^  <  n^ ,  j»l,2,...,k. 

Ve  now  establish  the  following  result,  which  may  be  inter¬ 
preted  as  a  duality  theorem.  This  will  prove  useful  in  some  of 
the  subsequent  material. 

Theorem  2.3.  If  f*(x;a)  «  a,  0<a<l,  has  at  least  one  solution  in 
a,  then 


*«ci)  -  ln£ 


ja|f*(x;a)  ■  ctj. 


(2.12) 


If  f*(x;a)  >  a  for  all  a,  then  a^.^  •  0. 

Proof.  Let 

c  ■  inf|a|f*(x;a)  >  <*j  •  (2.15) 

The  infimum  in  (2.15)  is  attained.  Thus,  there  exists  a  pQ  such 


t 


that  c  •  h(pQ).  If  f*(xja)  >  a  for  all  a,  let  pA  -*■  0,  i-1,2,.... 
Then  h(p)  0,  since  h(0)  *  0  and  h(p)  is  continuous,  and 

** tx)  • 

•  -* 

Now  assune  there  is  at  least  one  a  with  f  (x;a)  -  a.  Then 
f*(x;ag£-j)  >  a  and  therefore  c  <  •  If  c  <  ag(£j»  then 

c  ■  hCpQ)  and  f*(x;c)  ■  a,  which  is  a  contradiction. 

Renarks .  Again,  only  the  continuity  of  h(p)  was  used  in  the 
proof  of  Theorem  2.3.  Under  the  hypotheses  of  the  Corollary  to 
Theorem  2.2,  for  a  series  system,  agf£)  is  the  solution  in  a  of 

f* (x; a)  -  a  .  (2.14) 

The  general  theory  described  in  this  section  applies  as  well 
to  what  is  known  as  systems  with  repeated  components  (see,  e.g., 
Harris  and  Sons  C1973)).  For  such  systems,  there  are  1  <  m  <  k 
unknown  parameters  P i »P2  * ' *  * since  the  "repeated  components” 
are  assumed  to  have  identical  failure  probabilities.  This 
assumption  permits  the  experimenter  to  regard  the  data  as 
^ni*Yi^*  i*1*2* •••»■»  a»d  employ  the  previous  results. 

For  example,  if  a  series  system  of  k  components  has  otj  of 
one  type,  <*2  of  a  second,  ...,  oB  of  an  mth  type,  then 

«1  o-  o,  k 

Mp)  ■  Px  P2  •••  PB  *  “i  ’  ki  • 

3.  Counterexamples 

In  this  section  we  restrict  attention  to  series  systems  and 
employ  the  ordering  function 

k 

g(x)  ■  n  (n. -x.)/n.  , 
i-1  x  1  1 
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introduced  following  Theorem  2.1. 
cese  the  reliability  function  h(p) 
ization  we  have  for  (2.4) 


As  noted  previously,  in  this 
k 

■  Up,.  With  this  special- 
i-1  1 


tj  ■  nx (1-r^)  (S.l) 

and  for  each  fixed  0  <  ij  <  tj ,  0  <  ij  <  ,  . . . ,  0  <  i  ^  j  <  t  ^  j , 

r^1 

tj  -  n^U-r^/C  n  (hg-igD/ngl) »  2  <  j  <  k  .  (S . 2) 


whenever  g(x)  »  r  ,  1  <  n  <  s.  If  ■  ■  s,  then  r  »  and  a  «  0. 

1  SO 

For  ic  >  0,  X  >  0,  let 

Ip(*.X)  -  J*  tie"1(l-t)X“1dt  .  0  <  p  '  ,  (3.3) 

the  incomplete  beta  function. 

It  is  well-known  that  if  t  is  an  integer,  t  <  n,  we  have 

I  (?)  P"“V  -  I_ (n-t ,t+l)  .  (3.4) 

i-0  1  p 

In  Sudakov  (1974)  the  following  inequality  was  published. 

p~{* CX)  >  r^J  <  I  k  (nj-tj.tj+l)  .  (3.5) 

11  pi 
i-1  x 

This  inequality  and  generalizations  of  it  were  further 
studied  in  Harris  and  Sons  (1980,1981).  (3.5)  implies 


f* (*;*)  <  1>(nl-t1,t1*l)  , 


hence  its  usefulness.  However,  as  we  now  establish,  (3.5)  is  not 
universally  valid,  as  was  claimed  in  Sudakov  (1974). 

Let  (XpXj)  ■  (Xj.O)  and  let  (n1(n2)  -  (n^nj).  Then 
f(x)  •  (nj-Xjj/nj  and  tj  -  Xj.  Consider  P-{g(X)  >  r^j.  If 


,M'  T  -  T|r~' 


I  p. 
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•  (1»«),  0  <  *  <  1,  we  have 

Pp{*CX>  i  r«}  ■  P.{t"2-X2)/,,2  2 r.}  • 

since  P-^-oj  ■  1.  by  C2-D-  Consequently, 

p;{*w)  »  '»}  -  p.{x2  i  V1-1.’} 

■  P.{X2  :  J"l<I-rn)}  • 

Since  «  (n^-x^/n^ 

Pp{«C»  i  *.}  ■  P.{X2  i  2*l}  • 

Thus  fron  (3.4), 

Pp{*CX)  >  rB]  -  laC2Cn1-x1).  2x^1)  . 

The  Sudakov  inequality  implies  that 

Ia(2(ni-Xi),  2x^*1)  <  Ia(n1-x1,  Xj+1) 

or 

Ia(2nlra*  ^iCi-V1)  1  VnlV  niC1'rm)+1)  *  (3'6) 

Let  h2(t;n2,rn)  and  hjCtjn^r^)  denote  the  beta  density  functions 
corresponding  to  the  left  and  right  hand  side  of  (3.6),  respec¬ 
tively.  Then,  provided  n^r^  >  1.  there  is  an  e  >  0  such  that 

h2(t;n2,ra)  <  hjCtjnj.r^)  0  <  t  <  e,  1-e  <  t  <  1  . 

This  implies  that  h1(t;n2,rjB)  and  h2(t;n2,rn)  intersect  in  at 
least  two  points.  If  t*  is  such  an  intersection,  setting 
hj  (t;nj ,rB)/h2(t;n2,rB)  •  1  gives 
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n.r  n.fl-r  ) 

t  1  *(l-t)  1  * 


e(nl»ra)  >  0  . 


Thus,  for  1  <  a  <  s,  there  ere  exactly  two  such  intersections. 

Therefore  there  is  a  such  that 

o 


li  (niT*>  n! C1_rm) ^l)  -  I2  Cn2rm,  , 

o  o 

for  z  >  zQ, 

IzCnlr«’  <  I9^n2Tm$  n2(1-r*)  +  1) 


and  for  z  <  z  . 

o 

IZ(nlVnlC1'r«)*1)  >  1z(n2ru*  n2(1’rm)+1)  * 

Thus  for  z  >  zQ,  (3.6)  is  violated.  (3.6)  was  used  as  a  lenaa 
by  Sudakov  (1974)  to  prove  the  inequality  (3.5).  This  lemma  was 
also  employed  in  Harris  and  Soas  (1980,  1981).  It  is  the  falsity 
of  this  leaaa  which  invalidates  (3.5). 

Table  1  provides  some  illustrations  of  the  violation  of 
(3.5)  for  k  ■  2  and  selected  values  of  (n^,n2),  (Xj,x2).  The 
smallest  value  of  P^P2  for  which  this  violation  occurs  is  also 
given  in  the  table,  where  it  is  denoted  by  a*.  In  addition, 
f*(x;a*)  is  tabulated.  Thus  for  a  <  f*(x;a*),  (3.5)  is  valid. 

The  calculations  were  made  by  means  of  a  FORTRAN  program. 

Note  that  for  (nj,n2)  ■  (5,5)  and  (Xj,x2)  ■  (1,1),  the 
inequality  was  not  violated. 
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Table  1.  The  Smallest  a,  a*,  and  f*(x;a*) 


Cftj.nj) 

(X  j i Xj) 

a* 

f*(x;a*) 

(5,5) 

(1.1) 

1.0000 

1.0000 

(5,5) 

(5,3) 

.7454 

.9998 

(S,10) 

•  d.O) 

.8798 

.8909 

(5,15) 

(0,3) 

.8698 

.8791 

(5.30) 

d.O) 

.8498 

.8467 

4.  The  Theory  of  Key  Test  Results 

If  for  nj  <  n2  <  ...  <  n^,  (x^ , x2 , . . . » xfc)  *  (Xj ,0, . . . ,0) , 
k  >  2,  then  x  is  called  a  key  test  result.  Winterbottom  (1974) 
asserted  that  subject  to  x^  <  f(k,nj),  where  f(k,nj)  is  the 
solution  in  f  of 

njk-f-1  -  ktCn1-f)n5”1]1/k  ,  (4.1) 


we  have  a  is  the  solution  in  a  of 
8(x) 


Ia^nl"*l’ 


0  <  a  <  1 


(4.2) 


This  would  imply  the  inequality  (3.5),  which  we  have  disproved 
in  Section  3. 

As  we  subsequently  establish,  the  error  in  Winterbottom ' s 
(1974)  result  is  a  consequence  of  falsely  concluding  that 
f(k,nj)  depends  only  on  n^.  It  is  easy  to  be  led  to  this  conclu¬ 
sion  on  intuitive  grounds,  since  (n^-x^,  n^,...,  n^)  would  seen 
to  be  a  less  favorable  experimental  result  than  (n^-x^ ,n2> . . .  ,nk) , 
whenever  n^  >  n^  for  at  least  one  index  i,  2  <  i  <  k.  However,  we 


will  now  establish  that  a  "modified"  key  test  result  holds  for 
Xj  <  f (k.nj.nj, . . .,nk) ,  where  f Ck .nj.nj, . . . ,nk)  is  the  solution 


in  f  of 


i,  -f-1*  l  n.  •  kC(n.-f)  1!  n,]1/k  . 

1  1-2  1  1  i-2  x 


(4.5) 


Theorea  4.1.  If  nx  <  n2  <  . . .  <  nk  and  x  »  (Xj ,0, . . . ,0) ,  with 
X|  ^  f (k j n ^ t a 2 , • • • i ^  n^,  then 

(4.4) 

Jk  &  X 

'i 


P-{fCX)  >  ri 


}  <  1  k  (•!-*!.*!♦*>  » 
n  P< 

i-1 


where  g Cxj.O , . . . ,0)  -  r ^ . 
k 

Proof.  For  7  x .  <  n, ,  from  Marshall  and  Olkln  (1979,  p.  78), 

Tt -  i-1  1  X 

H  (u.-x.)  is  a  strictly  Sehur-concave  function  of 
i-1  11  k  k 

(n1-x1,n2~x2, . . . ,nk-xfc) .  Thus,  if  (nj-Xj)  is  fl*ed»  ^ni”xi 


Is  minimized 


ft  — 

at  (nj  -  x±,n2,.. .,nk) .  Equivalently,  ^  (n±-x1) 


is  minimized  for  vectors  of  the  type  x  -  (  J  x. ,0,...,0) 
k  ’  -  -  i-1  1 

J  x.  is  fixed. 

t _ a  * 


when 


i-1 


Let  z  -  (tj ,x2»  ♦  •  • » *j[)  with 


\  x.  <  *i  <  «i- 
i-1  x  1  x 


Then 


ft  ft 

I  (n.-x  )  >  n.-x.  ♦  I  n£  . 
i-1  11X1  i-2  1 


(4.5) 


For  each  fixed  value  of  J  (n.-x.),  we  have 

i-1  x  1 


ft  ft  ft  ft 

n  (n. -x. )  >  (n.  -  l  x.)  n  n.  >  (n.-x.)  JI  n.  .  (4.6) 
i-1  1  1  ‘  1  i-i  x  i-2  1  x  x  i-2 


In  order  that 
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{J|  -  j/"!-*!')  '  {Il1!1t"l-*i)  i  j,1"!-!5}-  (4’7) 


k 

we  Bust  heve  x.  ■  x,  *  . »  x.  ■  0.  Mote  that  If  J  *1  —  nl* 

*  i-1 

*  «  n  v  <  j.,.,.,*,  <  n,  ,  the  two  eets  cannot  coincide,  because 

1  1*  2  2  k  k 

(0,n2>< • • ,1^)  is  in  the  right  hand  set  but  not  the  left.  Fro« 

(4.6)  it  follows  that 


{*| -  “l'vj/i} e  ^  tni'xi:i (4,) 

k  k 

Equality  holds  if  max  II  (n^t^  *  (n1-*l)  n  ni  when  It"!"*!5  * 
k  i-1  1*2 

(n.-x  )  ♦  l  n.-l.  Pron  the  arithaetic-geonetric  Bean  inequality 
11  i-2  1 

this  is  true  whenever 


k 

n.-x.*  I  n.-l 
1  1  i-2  x 


(hj-Xi) 


k 

n  n. 
i-2  1 


(4.9) 


Note  that  equality  in  (4.7)  nay  still  hold  if  (4.9)  is 

k 

violated  since  (n.-x.+  J  n,-l)/k  nay  not  be  an  integer  or  nay  be 

1  1  i-2  x 

bigger  than  sobs  of  the  n^,  i-l,2,...,k.  Thus  if  Xj  is  the 
saallest  Xj  value  for  which  equality  holds  in  (4.9),  then 
f (k,n^,n2, . . .  >i»k)  *  ^1* 

If  Xj  <  f(k,nlfn2, ...,nk),  then,  fron  the  above. 


f*(x;a)  -  sup  P<  J  Y,  >  n.-x.  ♦  J  n.>  .  (4.10) 

k  'i-1  x  "  i-2  1 

H  p.-a 
i-1  1 

Writing  (4.10)  as  an  iterated  sub  and  noting  that 
It(n-x,x*l)  is  a  decreasing  function  of  n  for  fixed  x,  we  have 
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where  the  are  Independent  Mnoaial  random  variables  with  para¬ 
meters  Cnj.Pj^).  i-2,...,k.  Writing 


k  k  n. 

1  i-2  1  1.1  J.l  *> 


where  the  Y^  are  independent  Bernoulli  random  variables  with 

parameter  p. ,  a  result  o£  Pledger  and  Proschan  (1971)  may  be 
1  k  n 

employed  to  show  that  the  upper  tail  of  £  £  Y,.  is  a  Schur- 

i-1  j-1  13 

convex  function  of  (-lnpj,  -lnp^,...,  -lnpj,  -lnp2, . . ., -lnp2, . . . , 
-lnp^, . . .,-lnpk)  and  therefore  f*(x;a)  ■  I^Cnj-Xj,  Xj+1),  as 
required . 

As  is  discussed  below,  (4.3)  may  have  no  solutions.  In  such 
cases,  and  in  general,  it  is  possible  to  strengthen  (4.3). 


Corollary.  For  each  f,  form  the  vector  t  «  (*l,*2, . . . ,*k)  from 

n  ■  Cnj,n2, . . . ,nk)  by  continually  reducing  the  maximum  (s)  until 

the  subtractions  total  f*l,  f  >  0.  Denote  by  f  •  (k,nj ,n2,  . . . ,nk) 

the  first  f  for  which 

k  k 

H  x.  >  (n.-f)  H  n  . 
l-l  1  ’  1  i-2  1 


Then  (4.4)  holds  for  x^  <  f '(k.n^ ,n2» . • .»nk) . 

Proof.  The  proof  proceeds  exactly  as  for  Theorem  4.1' by  notins 

k  k 

that  I  maximixes  n  r,  subject  to  0  <  r,  <  n.  and  T  r.  - 
k  i-1  1  1  “  1  i-1  1 

T  n.-f-l.  This  follows  since  z  is  majorixed  by  r  and  the 
i-1  1 

product  is  strictly  Schur-concave . 


Remarks .  For  0  <  f  <  n.. ,  the  right  hand  side  of  (4.3)  is  concave 
decreasing.  The  left  hand  side  exceeds  the  right  hand  side  when 
f  *  nj .  If  the  left  hand  side  is  less  than  the  right  at  f  ■  0, 
there  is  exactly  one  solution  f,  0  <  f  <  n^ .  If  not,  there  are 
no  solutions.  There  is  always  a  solution  if  n^  *  n2  ■  ...  »  nk> 

Proa  the  Corollary  following  Theorem  4.1,  x1  *  0  satisfies  (4.4). 

If  •  n2  «  ...  *  nk,  (4.3)  reduces  to  (4.1)  which  is 
Ninterbottoa 's  (1974)  condition.  However,  s  should  be  replaced 
by  s+1  in  his  fornula,  which  also  has  a  sign  error.  As  an 
example,  for  k  ■  2,  n^  ■  n2  ■  SO,  from  Winterbottoa  (1974),  (4.4) 
is  stated  to  hold  for  xl  <  17  or  Hj-Xj  >  33.  However,  33 «S0  <  41*41, 
and  therefore  (4.4)  only  holds  for  x^  <  13  or  n^-x^  >  37,  as  the 
Corollary  to  Theorem  4.1  shows,  or  the  solution  of  (4.3),  which 
gives  f (2 , 50 , 50)  -  13.14. 

The  dependence  of  f  on  n  may  be  seen  by  considering  an 
example.  Let  k  »  2 ,  nj^  »  5,  n2  ■  10 .  Then  from  the  Corollary 
following  Theorem  4.1,  (4.4)  only  holds  for  Xj  ■  0,  whereas  for 
Uj  *  n2  «  5,  it  holds  for  x^  »  0,1,2,  and  3.  Thus  the  case  of 
equal  n^,  i«l, 2,..., k,  does  not  give  the  minimal  f.  In  fact,  it 
may  be  seen  that  if  n^  >  2n^,  then  (4.4)  holds  only  for  x^  »  0. 

5.  Concluding  Remarks 

From  Table  1,  it  seems  reasonable  to  conjecture  that  (3.S) 
is  valid  for  those  values  of  a,k,it  likely  to  arise  in  practice. 

The  authors  are  continuing  to  investigate  the  problem  and  hope 
to  report  more  precise  conditions  for  the  validity  of  (3.5)  in 
subsequent  work. 
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